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^ Link weights are the principal parameters of shortest path routing protocols, the most com- 

monly used protocols for IP networks. The problem of optimally setting link weights for 
unique shortest path routing is addressed. Due to the complexity of the constraints involved, 
there exist challenges to formulate the problem properly, so that a solution algorithm may 
be developed which could prove to be more efficient than those already in existence. In 
this paper, a novel complete formulation with a polynomial number of constraints is first 
introduced and then mathematically proved to be correct. It is further illustrated that the 
formulation has advantages over a prior one in terms of both constraint structure and model 
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size for a proposed decomposition method to solve the problem. 
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Shortest path routing protocols such as OSPF (Moy, 1998) are the most widely deployed and 

commonly used protocols for IP networks. In shortest path routing, each link is assigned a 

weight and traffic demands are routed along the shortest paths with respect to link weights, 



based on a shortest path first algorithm (Bertsekas and Gallager, 1992). Link weights are 
hence the principal parameters and an essential problem is to find an appropriate weight set 
for shortest path routing. 

A simple approach to set link weights is the hop-count method, assigning the weight of 
each link to one. The length of a path is therefore the number of hops. Another default 
approach recommended by Cisco is the inv-cap method, setting the weight of a link inversely 
proportional to its capacity, without taking traffic conditions into consideration. More gen- 
erally, the weight of a link may depend on and be related to its transmission capacity and 



1 Some results in this paper have been presented on the 4th International Conference on Networking 
(ICN 2005), Reunion Island, April 2005, INFORMS Annual Meeting 2006, Pittsburgh, PA, USA, November 
2006, and the 6th International Congress on Industrial and Applied Mathematics (ICIAM 2007), Zurich, 
Switzerland, July, 2007, respectively. 



its traffic load. Accordingly, a problem worth investigating is the task of finding an optimal 
weight set for shortest path routing, given a network topology, a projected traffic matrix 
(Feldmann et al., 2001), and an objective function. 

The problem has two instances, depending on whether multiple shortest paths or only 
a unique routing path from an origin node to a destination node is allowed. For the first 
instance, a number of heuristic methods have been introduced, based on genetic algorithms 
(Ericsson et al., 2002) and local search methods (Fortz and Thorup, 2000). For the second 
instance, the Lagrangian relaxation method (Lin and Wang, 1993) and a local search method 
(Ramakrishnan and Rodrigues, 2001) have been proposed. These methods have also resulted 
in good routing performances by being tested in particular data sets. On the other hand, 
with these heuristic methods, the problem is not completely formulated and so in general 
is not optimally solved. It would be desirable if in average cases optimal solutions could be 
obtained for reasonably large data instances from real-world applications. 

From a management point of view, unique-path routing uses much simpler routing mecha- 
nisms and allows for easier monitoring of traffic flows (Ameur and Gourdin, 2003). Therefore, 
this paper focuses on the unique-path instance. The problem is referred to as the unique 
shortest path routing problem. It is a reduction of the integer multi-commodity flow problem 
(Ahuja et al., 1993). 

The problem has been well studied and efforts have been made to formulate the problem 
mathematically (Ameur and Gourdin, 2003; Zhang and Rodosek, 2005). Mathematical mod- 
els have also been developed for other related problems such as network design and routing 
problems (Bley and Koch, 2002; Holmberg and Yuan, 2004). Most of these models have 
formulated the corresponding problems completely, whereas they have been either demand- 
based or path-based and have left space to further explore the structure properties of the 
problems, from which more efficient solution methods may be derived. 

The main goal of this paper is to mathematically model the problem, which would yield a 
new exact solution approach for real-world applications in average data instances. In Section 
2, the problem definition is first specified. Two different complete formulations, a new one 
and an existing one for comparison, are introduced in Section 3. The new formulation is then 
mathematically proved to be correct in Section 4. Differences between the two formulations 
on both constraint structure and model size are discussed in Section 5. Conclusions are 
drawn in Section 6. 



2. Problem Definition 

The unique shortest path routing problem is specified as follows: 
Given 

• A network topology, which is a directed graph structure Q = (A/", C), where 

— M is a finite set of nodes, each of which represents a router; and 

— £ is a set of directed links, each of which corresponds to a transmission link; (For 
each (i,j) G C, i is the starting node, j is the ending node, and c^- > is the link 
capacity.) 

• A traffic matrix, which is a set of demands T>; (It is assumed that there is at most one 
demand between each origin-destination pair. For each demand k G T), Sk G H is the 
origin node, tk G H is the destination node, and dk > is the required bandwidth. 
Accordingly, S is the set of all origin nodes, T s is the set of all destination nodes of 
demands originating from node s G S, and T> s is the set of all demands originating 
from node s G S.) 

• Lower and upper bounds of link weights, which are positive real numbers w m i n and 
to mas , respectively; and 

• An objective function, specifically, to maximize the sum of residual capacities, 
Find an optimal weight set -uv,-, (i,j) G £, subject to 

• Flow conservation constraints: for each demand, at each node, the sum of all incoming 
flows (including the demand bandwidth at the origin node) is equal to the sum of all 
outgoing flows (including the demand bandwidth at the destination node); 

• Link capacity constraints: for each link, the load of traffic flows transiting the link does 
not exceed the capacity of the link; 

• Path uniqueness constraints: each demand has a unique routing path; and 

• Path length constraints: for each demand, the length of each path assigned to route 
the demand is strictly less than that of any other possible and unassigned path to route 
the demand. 
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Figure 1: Illustration of the path length constraints 

3. Problem Formulation 

In this section, the problem is mathematically formulated from two different perspectives, 
based on the study of the problem properties. For comparison, an existing model (Zhang 
and Rodosek, 2005) is introduced in more detail first. 

3.1 A Demand-Based Model 

According to the characteristics of unique shortest path routing, the routing path of a demand 
is the shortest among all possible paths. For each link, the routing path of a demand either 
traverses the link or not. 

Figure [1] illustrates the relationships between the lengths of the shortest paths and link 
weights. Paths in thick lines are routing paths. The path length and path uniqueness 
constraints require that the length of the unique shortest path to route a demand is less 
than that of any other possible path from the origin node to the destination node. 

With respect to the constraints, there are three scenarios to be considered: 

• If the routing path of demand k traverses link (i,j), the length of the shortest path 
from node Sk to node j is the length of the shortest path from node Sk to node i plus 
the weight of link (i,j); 

4 



• If the routing path of demand k does not traverse link (i,j) but visits node j, the 
length of the shortest path from node s k to node j is strictly less than the sum of 
the length of the shortest path from node s k to node i and the weight of link (i,j)', 
(Otherwise, there would be at least two shortest paths to route demand k.) 

• If the routing path of demand k neither traverses link (i,j) nor visits node j, the length 
of the shortest path from node s k to node j is less than or equal to the sum of the 
length of the shortest path from node s k to node i and the weight of link (i,j). 

Based on the above observations on the relationships between the length of a shortest 
path and the weights of links that it traverses, the problem can be mathematically formulated 
as a demand-based model (DBM) as follows (by defining one routing decision variable for 
each link-demand pair): 

Routing decision variables: 

x^e{0,l},\fkeV,y(i,j)eC (l) 

is equal to 1 if and only if the routing path of demand k traverses link (i,j). The number of 
this set of variables is |Z>||£|. 
Link weight variables: 

Wij e [w min ,w max ],\/(i,j) e C (2) 

represents the routing cost of link (i,j). The number of this set of variables is \C\. 
Path length variables: 

It g [0,+oo),VseS,VieA/' (3) 

denotes the length of the shortest path from origin node s to node i. Obviously, if* is the 
length of the shortest path to route demand k G T> and l\ = 0, Vs G S. The number of this 
set of variables is |«S||A/"|. 

Flow conservation constraints: 

■1, if i = s k 
1, if i = t k 

h-.(h,i)ec j-(i,j)ec I 0, otherwise 

The number of this set of constraints is |X>| |A/"| . 
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J2 A- E 4 = < !> ** = ** ykeV^ieN (4) 



Link capacity constraints: 



J2 dkX iJ ^ c ii> V (^i) e£ 



(5) 
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The number of this set of constraints is \C\. 

Path uniqueness constraints: under the combined restrictions of the flow conservation 
constraints and the path length constraints, the constraints are satisfied automatically. 

Path length constraints: 



If < If + Wij 







,VkeV,V{i,j)eC (6) 

If > If + Wij - M (l - 4) 

where e and M are appropriate constants with < e <C M. The number of this set of 
constraints is 2|Z>||£|. By enumerating all possible values of the routing decision variables 
x k ^k G V, (i,j) G £, it can be verified that the linearized constraints (jSJ) are identical to 
those constraints (J7j) originally presented in logic forms, as illustrated in Figure [TJ 
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Objective function: 



which is equivalent to 
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(i,j)eC kev 
Accordingly, the resulting complete model is presented as 

Optimize flSJ) 



DBM: 



Subject to (@D, ©, ©, ©, ©, © 



In the following, an equivalent model to DBM, which will be used to verify the correctness 
of DBM in Section HJ is derived. 



A necessary condition of the unique shortest path routing problem is the sub-path opti- 
mality requirement, which states that any sub-path of a routing path is still a unique shortest 
path (Bley and Koch, 2002). Specifically, given an origin node s G S and a node i G M 
where i ^ s, it requires that all demands originating from s and visiting i use the same 
incoming link to i. Mathematically, it is formulated as 

V msxx k hi <l,i^s,\/seS,\/ieN. (9) 

h:(h,i)eC 

The above constraints can be alternatively presented with linear constraints, by intro- 
ducing a new variable yfj,s G S,(i,j) G C: 

^•>4,VfcGP s ,V(z,j)e£ and J2 y s hi <l,i^s,\/seS,\/ieN. 

h:(h,i)e£ 

In DBM, the sub-path optimality constraints are not explicitly included. In the following, 
it is proved that the constraints are implied by the path length constraints. 

Proposition 1 The path length constraints in DBM imply the sub-path optimality con- 
straints. 

Proof. Suppose there are two demands fei, ki G V, and s/ Cl = Sk 2 = s. Assume they would 
use two disjoint paths to traverse from node u to node v. Demand k\ would use path 

Pj = (J1J2) ->■ (J2J3) ->...-> (jn-l,jn), (ji,jl+i) e C, I = 1, . . . , n - 1, 
where j\ = u and j n = v, and demand &2 would use path 

Pi = {hi I2) ->■ (*2, h) ->■■•• -> (im-l, »m), (*g, VfO E £,q = 1, . . . ,171 - 1, 

where %\ — u and i m = v. 

Then, by the definition of the routing decision variables, 
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fci _ 1 „„^l „te 



1 and xf in = 0, V(ji, ji+i) ePj,l = l,...,n-l. 



JUl+i JUl+i 



As a result, according to constraints (|7j), on one hand, since VA; G V,W(i,j) E C, x 
> l^ = ll k + u>„, considering demand k\, 



ll + l? 

u Jn 




l S u + li^+W jn _ ljn 




l u + ll- 2 + W jn _ 2jn . 


1 ' Jn — ljn 


l 8 u + l\ + w jlj2 + • • 


1 Jn — ljn 


l S u + l Pj - 





(10) 



On the other hand, since Vfc G V,V(i,j) G C, x^ = A J2h-thj)ec x hj = 1 =^ If < ^t k + u!ij 
and Xy = =>- /^ fc < Z* fe + Wy, considering demand &2, 

I s = I s + H 
'u '« ^ 'in 

< ^ + ^„_ 1 +^in- 1 in 

< ^ + ^„_ 2 + ^in_2in-l + ^in-lin 

< Z* H- J£ + w jin + ■■■ + w jn _ ljn 

= l S u + lPr (11) 

Contradiction then follows between ffTUl) and (fTTj) . which means the two demands cannot 
be routed over two different paths between two shared nodes. It is hence proved that the 
sub-path optimality constraints are satisfied. D 

By Proposition [I] DBM is equivalent to the following model DBM'. 

_ Q *T)t%m% z C (181) 

DBM ' : Subject to ^MMMMMM 

3.2 An Origin-Based Model 



In Section 13.11 the unique shortest path routing problem is formulated as a demand-based 
model, which defines one routing decision variable for each link-demand pair. Based on the 
study of solution properties of the problem, it can be found that all routing paths of demands 
originating from the same node constitute a tree, rooted at the origin node. Accordingly, a 
more natural formulation for the problem is to define one routing decision variable for each 
link-origin pair. For example, in Figure [2], instead of defining three routing decision variables 
for link (i,j), one for each of the three demands sharing the same origin node s, the new 
formulation defines only one routing decision variable for link (i,j), associated with origin 
node s. 

Based on the above observations, an origin-based model (OBM) for the problem is for- 
mulated as follows. 

Routing decision variables: 

^.e{0,l},VaeS,V(i,j)e£ (12) 

is equal to 1 if and only if the routing paths of at least one of the demands originating from 
node s traverse link (i,j). The number of this set of variables is |«S||£|. 
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Figure 2: Illustration of the origin-based model 



Auxiliary flow variables: 



/?.e[0,+oo),V 5 eS,V(i,.7)e£ 



(13) 



represents the load of traffic flows originating from node s and traversing link (i,j). The 
number of this set of variables is |<S||£|. 
Link weight variables: 

Wij € [ty m in,Wmax],V(i, j) G C (14) 

represents the routing cost of link (i,j). The number of this set of variables is \C\. 
Path length variables: 

l*e [0,+oo)yseS,VieAf (15) 

denotes the length of the shortest path from origin node s to node i. In particular, l s s = 
0, Vs G S. The number of this set of variables is |«S||A/*|. 

Flow conservation constraints: for each tree, at the root node, the difference between the 
sum of outgoing flows and the sum of incoming flows is the sum of bandwidths of all demands 
originating from the root node; at the destination node of each demand originating from the 
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root node, the difference between the sum of incoming flows and the sum of outgoing flows 
is the bandwidth of the demand; and the sum of incoming flows is equal to that of outgoing 
flows at any other node. 

{— d s , if i = s 
d k , iU = t k ,VkeV s ,VseS,VieM (16) 

0, otherwise 



where d s = X]fcez> ^fc- The number of this set of constraints is |«S| 

Flow bound constraints: for each tree, the total flow load over each link does not exceed 
the sum of all demand bandwidths originating from the root node and it is equal to zero if 
no demand originating from the root node is routed over the link. 

%<ri s Y,d k ,V8eSMiJ)e£ (17) 

kev s 

The number of this set of constraints is |<S||£|. 
Link capacity constraints: 

£^< % ,V(i,j)e£ (18) 

seS 

The number of this set of constraints is \C\. 

Path uniqueness constraints: for each tree, the number of incoming links with non-zero 
flows is equal to zero at the origin node; the number of incoming links with non-zero flows 
is equal to one at the destination node of each demand originating from the root node; and 
the number of incoming links with non-zero flows does not exceed one at any other node. 



(19) 



The number of this set of constraints is \S\ 

Path length constraints: for each tree, the length of the unique shortest path to route a 
demand originating from the root node is less than that of any other possible path from the 
origin node to the destination node. 
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},WseS,\/(hj)eC (20) 



The logic constraints (120 j) can be linearized as follows: 

i] < n +wij - e l j] ^. - ^. I 5j e £ 

\h:(h,3)eC J 

i) > n + w l3 - m(i - 4.) j 

where £ and M are appropriate constants with < e <C M. The number of this set of 
constraints is 2|«S||£|. 
Objective function: 




which is equivalent to 



min £ E^ ( 22 ) 

(ij')6£ seS 



Accordingly, the resulting complete model is presented as 
Optimize (1221) 

OBM: S'w&ject to ([ISD, (HZD, dUD, (USD, (EH), (H2D, (O, ([HD, (USD 
4. Model Verification 

The correctness of OBM is verified in this section. The verification is divided into two steps. 

First, DBM is verified to be a correct model of the unique shortest path routing problem. 

Second, OBM and DBM are mathematically proved to be equivalent concerning both 
the feasibility and the optimality of the problem, which implies that OBM is also a correct 
model of the problem. 

Since a relaxation of the unique shortest path routing problem has been thoroughly 
studied and a corresponding demand-based formulation for the relaxed problem has been 
recognized as being correct, the verification of DBM is built on the correctness of the for- 
mulation for the relaxed problem, the integer multi- commodity flow problem (Ahuja et al., 
1993) with sub-path optimality condition. Also, the equivalence between OBM and DBM is 
demonstrated based on the proof of the equivalence between two corresponding models of 
the relaxation, which are denoted as RDBM and ROBM, respectively. 

Rr) r, M . Optimize flSJ) 

Subject to ®, ©, ©, © 

Optimize 



robm: Subject to ([xbd, (HZD, tfHD, (USD, flED, flU 
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As can be noted, RDBM is actually a relaxation of DBM', which is equivalent to DBM 
by Proposition [TJ 

4.1 Correctness of DBM 

To verify the correctness of DBM, that of the equivalent model DBM' is proved. Apparently, 
the difference between DBM' and RDBM lies in the path length constraints ([6]) and the 
additional link weight variables ([2]) as well as the path length variables ([3]) . In order to verify 
that DBM' formulates the unique shortest path routing problem correctly, constraints (J6]) in 
DBM' are proved to represent correctly the additional path length constraints. Specifically, 
the following two statements are demonstrated to be correct. In DBM', the path length 
constraints ([6]), combined with the flow conservation constraints (j4]), guarantee that: 

1. The routing path of each demand is a shortest path; and 

2. The routing path of each demand is a unique path. 

The two statements are verified in Proposition [2] and Proposition [3j respectively. Propo- 
sition [2] is proved by demonstrating that the length of the routing path of each demand is 
less than or equal to that of any other possible path. The proof of Proposition [3] is built 
on, with two lemmas, the satisfaction of the single-path requirement by the relaxed problem 
RDBM, followed by the proof that the uniqueness requirement is satisfied by DBM'. As 
the original logic constraints ([7|) are identical to the linearized constraints fl6]), the following 
proof is based on the original constraints. 

Proposition 2 Each routing path resulting from the solution to DBM is a shortest path. 
Proof. Assume for demand k, 

Pj = (jlih) -> (h,js) ->-..-> (jn-l,jn), (jljl+l) e C,l = 1, . . . ,n - 1, 
where j\ = Sk and j n — tk, is the assigned routing path, and 

Pi = («i, z'a) ->■ (*2, h) -»• ■ ■ • -> (w-i, i m ), (V Vt-i) e £, g = 1, . . . , m - 1, 

where %\ = Sk and i m = tk, is any other possible and non- assigned path from s^ to £&. 

Then, by the definition of the routing decision variables, a;* • = 1, V(ji,ji+i) € Pj, 
I = 1, . . . , n - 1 and 3(i q , i q+x ) G P h x k iqiq+i = 0, q e {1, 2, . . . , m - 1}. 
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As a result, according to constraints fl7|), on one hand, since V(i,j) G C, x% = 1 => Z^* 



Z, Sfe + w 



y 3 

Z fc = 5* = ^-1 + %-U, = ^i + W ilj2 + ■ ■ ■ + Wjn-Un = * P,- • (23) 



On the other hand, since V(z, j) G £, a£ = =>- l Sk < l s t k + m 



'./• 



£* = 7 Sfc < 7 Sfc + w - • < 7 Sfc + Wii + . . . + w . . = /p.. (24) 

It thus follows that I p. < lp i from (|23|) and (|24p . Hence, the length of the routing path 
Pj is less than or equal to the length of any other possible path Pi to route demand k. It is 
therefore proved that the routing path Pj is the shortest one. D 

In order to prove that DBM' satisfies the path uniqueness constraints, the constraints 
are demonstrated to be satisfied by the relaxed problem RDBM first. 

Lemma 1 No optimal solution to RDBM contains flow loops. 



Proof. Suppose x*\j,k G V, (i,j) G C is an optimal solution to RDBM, then all con- 
straints PJ, (jHJ), and (EJ) are satisfied by x*^, k G V, (i,j) G £. Specifically, 

{-1, if i = s k 
1, M = t k ,VfcGP,V*GA/-, (25) 

0, otherwise 

E d k x*% < c ij} V(i } j) G C, (26) 



kev 



and 



V maxx* k hi <l,i^s,VseS,VieAf. (27) 

h:(h,i)e£ 

Assume there would exist k\ G V and a loop C : (ji,J2) —>••••—>■ (j n -i,3n) ~~ ^ 

,jn+l),{ 

Obviously, 



(jn,jn+i),(ji,ji+i) G £,i = l,...,n and j„+i = ji, such that £*^ +1 = l,i = l,...,n. 



Let 



E ^- E ^t^O^zGJj!,...,^}- (28) 

h-.(h,i)ec j-.(i,j)eC 



'■' I 0, otherwise , V/c G i^, V^, jj G >C. ^Jj 
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Then, y*% G {0, 1}, Vfc G V, V(i, j) G £ and y% +1 = 0, V&,*+i) E C,i = 1, . . . ,n. In 
addition, according to ( 125|) . fl28 |) . and f l29|) . 

{-1, iii = s k 
1, if z = t fc ,Vfc€P,V^GAT, 
0, otherwise 



E d *y* 



according to (1261) and ( l29i) . 

^2dkV*ij < E^ x *?? < c i:j ,V(i,j) G £, 

k&D keV 

and according to (127j) and (j2"5|) . 

E maxj/t, < E ™ XI 1. < M 7^ ^ 5 Vs G <S,Vi G W. 

h:(h,i)eC h:(h,i)eC 

Hence, y*^,k G T>, (z,j) G £ satisfies all constraints (TJJ, (EJ), and Q). It is therefore a 
feasible solution to RDBM. Furthermore, according to (I2§|) . since V/c eD,4> 0, 

E E^t- = E *<+ E^*§ 

(i,j)e£ kev (i,j)ec\cvkev-.k^ki (i,j)ec 

tk 

■i-j 

(i,j)eC\C\/keV:k^k t 

= / "fc^ ij 

(i,j)eC\C\/k£T>:kj£ki 

< Yl d k x*ij+ Y d h x *ij 

{i,j)ec\cvkeV:k^h (i,j)eC 

= E E 4x t- 

Then, £%„■, k G P, (i, j) G £ could not be an optimal solution to RDBM, which results in 
contradiction. It hence proves that no optimal solution to RDBM contains flow loops. □ 

Based on Lemma [U the path uniqueness constraints are proved to be satisfied by RDBM 
in Lemma |5J 

Lemma 2 The path uniqueness constraints are satisfied by RDBM. 

Proof. By Lemma [Q in RDBM, the flow conservation constraints fll]) at origin nodes are 
equivalent to: 

Y x hi = and Y x ij = M = s k , VA; G V. (30) 

h-.{h,i)ec j-(ij)e£ 
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Similarly, the flow conservation constraints fll]) at destination nodes are equivalent to: 

J2 a& = land J2 x k 3 = 0,i = t k ,Vk eV. (31) 

Since x k j G {0,1}, VA; G V,W(i,j) G £, constraints f|30|) restrict that there is one, and 
only one, outgoing link with a non-zero flow from the origin node of demand k. Similarly, 
constraints (13TJ) restrict that there is one, and only one, incoming link with a non-zero flow 
into the destination node of demand k. In addition, constraints fll]) guarantee that, at each 
intermediate node, the number of incoming links with non-zero flows is equal to the number 
of outgoing links with non-zero flows. Hence, for each demand k, the number of routing 
paths is no more than one. Therefore, the path uniqueness constraints are satisfied. □ 

Proposition 3 The path length constraints in DBM restrict that the resulting shortest path 
of each demand is a unique path. 

Proof. Since DBM' is a reduction of RDBM , the solution to the routing decision variables 
x k p keV, (i,j) G C of DBM' is also a solution to those of RDBM. 

By Lemma [2j there is only one routing path for each demand. Suppose for demand k, 

Pj = {h,h) -> (h,js) -►...-► (jn-i,jn), (jiji+i) e C, l = 1, • . . , n - 1, 

where j\ = Sk,j n = tk, is the assigned routing path, and 

Pi = {h, k) -> (l2, is) -> • • • -» («m-l, im), (V iq+i) E C, q = 1, . . . , m - 1, 

where %\ = Sk and i m = tk, is any other possible and non-assigned path from Sk to tk- 

Then, by the definition of the routing decision variables, X* ■ i = 1, V(j/, j/+i) G Pj,l = 

l,...,n-l. 

As a result, according to constraints (EJ), since V(i, j) G C, x k j — 1 =>• Z* fe = Z 4 Sfc + iu^-, 

£ = '2 = *£_i + w ^u = # + w nn +■■■ + u^U» = h r (32) 

As both Pj and Pj are paths between s^ and tk, they finally merge at one node. Assume 
it is node r and r = j p = i q ,p G {2, 3, ... ,n — l,n}, g G {2, 3, ... ,m — l,m}. Then, by the 
definition of the routing decision variables, x k r = 1 and x k „ = 0. 

In addition, on one hand, S/i:(/i,r)e£ x ftr < 1- On the other hand, Z]h:(h,r)e£ x L ^ x l_ir + 



X i q -ir — 1- Hence, Z^h:(h,r)€C X hr ~ *- 
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As a result, according to constraints (ED, 



'« m _l "■" W im-lim 

= K k + w riq+1 + ...+ w im _ lim 

< l i^ x + W i g -ir + Wri q+1 +■■■ + W im _ lim 

< 111 + W ni2 + ■■■+ W im _ lim 

= l Pi - (33) 

It follows that lp j < lp i from ( |32|) and ( 1331 . It is therefore proved that path Pj is the 
unique shortest path to route demand k. D 

Corollary 1 DBM is a correct model of the unique shortest path routing problem. 

Proof. By Proposition [2] and Proposition [3J DBM' is a correct model of the unique 
shortest path routing problem. Hence, as an equivalent model to DBM', DBM is also a 
correct model of the problem. □ 

4.2 Correctness of OBM 

By the proof of Lemma [2J the flow conservation constraints (jlj) are identical to: 

Y ^ = °> Y 4 = 1 > [ii = s k 

h-.(h,i)ec j-.(i,j)eC 

Y x hi = 1 > Y 4 = °> iU = f k 

h-.(h,i)ec _ j--(i,j)eC 

Y x m = Y x iP otherwise 
h:(h,i)eC j-.(i,j)eC 

Hence, RDBM is equivalent to the following model: 

Rr)RM /. Optimize flS} 

. Subject to (El , © , (ED , dU 

Lemma 3 In ROBM, the flow conservation constraints at origin nodes are equivalent to: 
Y ft = and Y H,= Y d ^ l = s ^ seS - 

h:(h,i)eC j'-(i,j)e£ k&V a 



> ,VkeV,\/ieN. (34) 
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Proof. On one hand, according to (I13p . if i = s, 

E /&>0,V*eS. 

h:(h,i)e£ 

On the other hand, according to (IT7|) and (fT^l) . if i = s, 



E /*< E (»*!>*) = £* E I^ = 0,VaG5. 

h:(h,i)eC h:(h,i)£C \ keV s J keV s h:(h,i)eC 

Hence, 

h:(h,i)eC 

It can then be derived directly from (fl6j) that 

E # = E^ = s ' VseiS - 

i:(ij)e£ A:eX> s 

The conclusion is therefore verified. D 

By Lemma [HI the flow conservation constraints ( fT6l) are identical to: 

E f^ = o, E /$ = <*•• iH = s 

h:(h,i)£C y.(ij)ec 

E /«- E # = *• if< = t*,V*eZ>. 

h-.(h,i)ec j--(i,j)ec 

E ^ ~ E $ = °> ° therwise 

h-.(h,i)ec j-(i,j)e£ 



>,VseS,VieAf, (35) 



where 4 = Efcex> a 4- 

As a result, ROBM is equivalent to the following model: 

RORM- °P timize d22D 

S'Mftjec* to (|35D, (HZD, (USD, (BSD. dI2D, (BSD 

In the following, ROBM and RDBM are proved to be equivalent concerning the feasibility 
of the relaxed problem. The proof is heavily based on the verification of the equivalence 
between the two respectively identical models, ROBM' and RDBM'. 

Proposition 4 There is a solution to ROBM if RDBM is feasible. 
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Proof. Suppose x*^, k G T>, (i,j) G C is a feasible solution to RDBM, and so to RDBM'. 
Then, all constraints f )34|) . (jSJ), and (jUJ) are satisfied by x*^-, k G V, (i,j) G £. 
Let 

y*£ = maxx*J,Vs G S,V(i,j) G £, (36) 

and 

/•£ = Y d k x%Vs G 5, V(i,j) G £. (37) 

fce£> s 

Obviously, y*ij G {0, 1} and /*£ G [0, +oo), Vs G <S, V(i, j) G £. 



According to (01, VA; G P, if z = s k , J2 h :( hi)e c x * k hi = ° and Y, r u mc x % = L Then ' 



according to (1371). Vs G 5, if i = s, 



Y f*M = Y Y dkX *hi = Y dk Y X *hi = °> 

h:(h,i)eC h:(h,i)eCkeT> 3 keT> a h:(h,i)eC 

and 

Y n* = Y Y dkX % =Y dk Y x *% = Y dk - 

r-(i,j)ec j:(i,j)€£k£T> a kev s j-.(i,j)ec kev a 

Hence, at original nodes, constraints (|35|) in ROBM' are satisfied by /**-, s G S, (i,j) G C. 

Also, according to (El, \fh G V, if i = t h , J2h:(h,i)&c x *hi ~ T,y.(i,j)ec x *ij = X and 
VA; G V Ski ,k ± k h E MMe ^l - Ej:(i, j} ec x *i = °- Then ; according to flST), Vfc, G 2>, if 
s = s kl and z = t kl , 



Y f*M~ Y f*ii = Y Y dkX *hi~ Y Y dkX * 



■fk 
ij 
h:(h,i)eC j-(i,j)eC h:(h,i)GCk£V s j:(i,j)££k£V s 



i^ i d k Y x * k hi ) - Y [ dk Y x % 

k&V s \ h:(h,i)eC I k£V 3 \ j:(i,j)eC 



Y dk [ Y x *hi - Y x % 

keV s \h:(h,i)eC j:(i,j)€C 



d h [ Y x * k hi~ Y 

.h-.{h,i)ec j-(i,j)ec 



x i 



+ E * E *u- E 

d, 






tfc r 



Hence, at destination nodes, constraints (1351) in ROBM' are satisfied by /*?-, s G <S, (i, j) G £. 

Similarly, according to (EH), VA; G £>, if i 7^ s fc and i ^ t k , Y.h:{h,i)&c x *hi-Y.j:(i,j)ec x * k ij = 

0. Then, according to ([37}, Vs G £>, V? G TV, if i ^ s and i $ V s 



'Sl 



E f*M- E f% = E J2 dkX * k ^- E E 4x t- 

h-.(h,i)ec j:(ij)e£ h:(h,i)eckev 3 j-.(i,j)eckev s 

= E U E *•!) - E (* E 

feel's \ h:(h,i)e£ I k£V 3 \ j:(i,j)eC 



™*fe 



y^ dfc y^ •'•'"/,; - > •'•' 



hi _ Z^ X ij 

keV s \h:(h,i)eC j:(i,j)GC 

0. 



Hence, at other nodes, constraints (135|) in ROBM' are satisfied by f**j, s G 5, (i, j) G £. 
According to (USD and ([37]), Vs G S,V{i,j) G £, 

Hi = E 4x t ^ E 4 K x t- = £^ x t- E 4 = ^*y E 4 - 

fcex> 3 kev s s " fce£> s fee£> s 

Constraints (PU) in ROBM' are then satisfied by $/*£, /*£, s G 5, (z, j) G £. 
According to (J3"7I) and (J5J), V(i,j) G £, 

En,- = EE d * x *S = E d * x t- ^ %• 

seS ses fcex> s fcex> 

Constraints flTHj) in ROBM' are thus satisfied by /*•.,•, s G 5, (i,j) G £. 

According to ([21]), VA; G £>, if i = s k , J2 h , {hti)£C x *hi = ° and so vfc e V i V(M) G £, if 
z = Sfe, a:**, = 0. Then, according to (13"oT) . Vs G 5, if i = s, 

E ^= E k x ***= E s = - 

h:(h,i)eC h:(h,i)£C " h:(h,i)eC 

Hence, at original nodes, constraints ([Til in ROBM' are satisfied by y**-, s G 5, (i, j) G £. 

Also, according to (1M|) . V/ej G X>, if i = t^,, J2h-(hi)ec x *hi = 1- Then, according to (15B|) . 
Vfc/ G X>, if s = sk and i = tk,, 



E*S \ V *fc v. \ V *fc; -1 

^^ = 2_, S x ^- 2_, x hi = l - 

h:(h,i)eC h:(h,i)€C h:(h,i)£C 



In addition, according to (jUJ), V^ E V, if i = t kl , Y,h:(h,i)ec max k&T> Sk x* hi < 1. Then, 
according to (1351) . Vfc; G X>, if s = s^, and i = t k 



1 1 



^ w = Z^ S x w < 1- 

h:(h,i)eC h:(h,i)£C 
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Thus, VA; G V, if s = Sk and i = tk, 

E y*^ = L 

h:{h,i)eC 

Hence, at destination nodes, constraints ( jTTJ|) in ROBM' are satisfied by y*?-, s G S, (i,j) G £. 
Similarly, according to 0, Vs G 5, Vi G A/", if i ^ s, Y2h-.(h,i)ec max fce£> 3 x *hi < 1- Then, 
according to fl36l) . Vs & S,Wi E Af, if z 7^ s and i ^ T s , 

E **« = E K**^ 1 - 

h:(h,i)eC h:(h,i)eC 

Hence, at other nodes, constraints f|T9|) in ROBM' are satisfied by ?/*?■, s £ S, (i,j) G £. 

Since all constraints flSSD, (H3), (USD, and (USD in ROBM' are satisfied by y*^ and /*£, s G 
5, (i, j) G £, it is a corresponding feasible solution to ROBM' and also to ROBM, of the 
feasible solution to RDBM, x* k ip keV, (i,j) G C. D 

Proposition 5 There is a solution to RDBM if ROBM is feasible. 

Proof Suppose y*ij,f*tj,s G S,(i,j) G C is a feasible solution to ROBM, and so to 
ROBM'. Then, all constraints (jMJ, (DID, (USD, and flT2} are satisfied by y%f%s G 
S,(i,j)eC. 

Let 



or 



/if 



/t 



y*% k ,VkeV,V(h,t k )e£. 



According to constraints f lT9|) at destination nodes, J2h-(ht )e£?/*/\t fc = IjVA; G XX Hence, 
VA; eV,3h: (h, t k ) G £, such that y* s h k t = 1. 



VA; G V,V(i,j) G £, if (z,j) 7^ (i,tk), x*^ is assigned as follows: 






Initialize x*J- <- 0, Vk eV,V(i,j) e £,(i,j) ^ (i,t k ) 
For fceD 

z G- 4 

Do find /i : (/i, z) G £, such that y*J* = 1 (38) 

Until i = Sk 

In the above assigning process (13"8"j) . at node i in each iteration of the inner loop, on one 
hand, according to (j3"o"|) . 

E /*J3> E n:,+4>4>o,vA;GX>, 

h-.(h,i)ec j-(t k ,j)e£ 
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and according to (TT7|) . 



Hence, 



E n?< E (v*s E <** = E ^ E ^.v^p- 

h:(h,i)eC h:(h,i)eC \ k'&V 8h J k '^ v s k h:(h,i)eC 



E y*S>0,V*e2>. 

h:(h,i)eC 

On the other hand, according to (!T9|) . 

h:(h,i)e£ 

Then, at node i in each iteration of the inner loop, 

E y'3 = i,V*e2>. 

h:(h,i)eC 

Therefore, at node i in each iteration of the inner loop, 3h : (h, i) G C, such that x* k hi 
V*hi = 1- Moreover, VA; G D, according to (I55|) . the process terminates at node Sk- 
According to (J551) . obviously, 

a:** G {0, 1} and *•* < y*J,V* G D,V(i,j) G £. 

Then, according to ffT9|) at original nodes, if i = s&, 



E ^ E v m Z = Wkel>. 

h:{h,i)£C h:(h,i)<=C 

In addition, the assigning process terminates at node Sk, Vfc G V, and so if i = S&, 

E s*J = l,VA;e2>. 

i:(*,i)e£ 

Hence, at original nodes, constraints (1511) in RDBM' are satisfied by x*y, fceD, (i, j) G £. 
According to (TT9l) at destination nodes, if i = t&, 

E *** = E V*JB = l,V*e2>. 

h:(h,i)eC h:(h,i)e£ 
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Furthermore, according to (1351) . if i = t k} x*f- = 0, \/k G V,V(i,j) G C. Then, if i = t 



i.i 



E x*^ = o,WkeV. 



Hence, at destination nodes, constraints ( |3lj) in RDBM' are satisfied by x*^, k G T>, (i,j) G C. 
According to (I38j) . at node i in each iteration of the inner loop, 

E *** = E **S = l,V*eX>, 

h-.(h,i)ec j:(i,j)e£ 

and at any other node i' G A/", i' ^ Sk,i' ^tk-, 

E *•«'= E 4 = °> vfc ^' 

h-.(h,i')ec r.(i',j)e£ 



Then, if z 7^ Sfc,2 7^ ^ 



E **&" E < = 0,V*e2>. 

h-.(h,i)ec j-{i,j)ec 

Hence, at other nodes, constraints fl34|) in RDBM' are satisfied by a;*y, k E V, (i,j) G £. 
According to (I35I) and (TT81) . 

E d fc x *£ = E E rffcX *t' - E-H? - c io^(hJ) e £. 

Hence, constraints (J3J) in RDBM' are satisfied by x*^-, fe G V, (i,j) G £. 

According to fl3_HD, a;**- < y*-. fc , Vfc G £>, V(i,j) G £. Then, according to ([T91, if i / s, 



E ^ XX *^- E £5*^ = E yli<l,V*G<S,ViGAA 

h:(h,i)£C h:(h,i)€C h:(h,i)eC 

Hence, constraints ([9]) in RDBM' are satisfied by x*^, fcGD, (i, j) G £. 

Since all constraints fl3jj, ©, and ([9]) in RDBM' are satisfied by x* k ij} k G V, (i,j) G £, 
it is a corresponding feasible solution to RDBM', and so to RDBM, of the feasible solution 
toROBM,y%f%seS,(i,j)eC. U 

Theorem 1 ROBM and RDBM are equivalent concerning the feasibility of the relaxed prob- 
lem. 

Proof. The conclusion is derived directly from Proposition H] and Proposition |5j D 
Based on the proof of the equivalence between ROBM and RDBM, the equivalence be- 
tween OBM and DBM, concerning the feasibility of the unique shortest path routing problem, 
is verified as follows. 
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Proposition 6 There is a corresponding solution satisfying the path length constraints in 
OBM, for each solution satisfying the path length constraints in DBM. 

Proof. Suppose x*^, k G T>, (i,j) G C, w*ij, (i,j) G C, and l* s t , s G S, i G M is a feasible 
solution to DBM and also to DBM' by Proposition [TJ Since RDBM is a relaxation of DBM', 
x*ij,k G V, (i,j) G C is then a feasible solution to RDBM and satisfies all corresponding 
constraints. 

Let 

y%. = maxx%V S GS,V(U)G£. 
Obviously, y* s {j G {0, 1}, Vs G S, V(z, j) G £. Also, Vs G S,\f(i,j) G £, there are three cases: 
• Case 1: *** = and Eh:(hj)eC^h j = 0, Vfc G P s ; 

1 =12-. ^ T> V^ 

:(fej)6£ A hj - x ! 



Case 2: x** = 0, VA; G P s and 3 A;, G P s , £ fc:(h ?)e £^fe = 1| 



• Case 3: 3ifej 6D S ,4 = 1. 

For Case 1, on one hand, Vi 6 P s , constraints (jBj) can be simplified as follows: 

If < Jf* + via and Zf > lt k + tw y - M. 

On the other hand, since y*?. = max fee -p s x*^- = and Yjh-(h j)ec y*M = ^> constraints (12~T1) 
can be simplified as follows: 

lj < It + Wij and lj > 11 + Wij - M. 

Hence, the simplified constraints of (12TT) are identical to those of ([6]) for Case 1. 
For Case 2, on one hand, constraints ([6]) can be simplified as follows: 

l] kl < l- kl + w^ - e and if 1 > Q l + Wij - M, 

k 



and V/c G V s such that ^2h-(hj)ec x *hj = 0> 



Zf < /" fc + w^ and Zf > /" fc + Wij - M. 

On the other hand, since y*|- = maxj. e p s x*^- = but 52ft-(h j)€cV*hj = 1> constraints (I2~T1) 
can be simplified as follows: 

lj < l! + w^ - e and lj > /• + w {j - M. 
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Hence, the simplified constraints of (|2"Tj) are identical to those of (E|) for Case 2. 

For Case 3, on one hand, according to the sub-path optimality constraints, Vfc £ P s , if 
0) Eh-(h j)ec x *hj = 0- Then, constraints (EJ) can be simplified as follows: 



X ij 



l/ 1 < l* 1 + Wij and I* 1 >l i h + Wij, 



and V/c £ £> s such that x*\j = 0, 

Zf < l- h + w^ and If > l- k + w {j - M. 

On the other hand, since y*^ = max^ e p s i*^ = 1 and Y^h(hj)ecy*hj = 1; constraints (I2TJ) 
can be simplified as follows: 

lj < It + w^ and I j > 11 + Wij. 

Hence, the simplified constraints of (|2"Tj) are identical to those of (E]) for Case 3. 

Since for all the three cases, y*tj,s £ S,(i,j) £ C results in the same path length 
constraints for OBM as those resulting from x*£.-, k £ X>, (z, j) £ £ for DBM, then, there is a 
corresponding feasible solution satisfying the path length constraints (12"T1) in OBM, provided 
that there is a feasible solution satisfying the path length constraints §§§ in DBM. D 

Proposition 7 There is a corresponding solution satisfying the path length constraints in 
DBM, for each solution satisfying the path length constraints in OBM. 

Proof. Suppose y*fj,f*tj,s £ S,(i,j) £ C, w*ij,(i,j) £ £, and l**,s £ S, i £ J\f is a 
feasible solution to OBM. Then since ROBM is a relaxation of OBM, y*^, f** j7 s £ S, (i,j) £ 
L is also a feasible solution to ROBM and so satisfies all corresponding constraints. 

According to the assigning process ( 13"5j) . let x*^, k £ V, (i,j) £ £ be the corresponding 
solution to RDBM. 

Vs £ S,W(i,j) £ £, there are three cases: 

• Case 1: y* s {j = and £ fc:(fc j )e c2/*jy = 0; 

• Case 2: y*^ = and Efc : (fcj) B c2/*&j = ^ 

• Case 3: ?/% = 1. 
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For Case 1, constraints (j2ip can be simplified as follows: 

lj < 11 + Wij and I] > 1° + Wij - M. 

According to ([38]), Vfc G V s ,x*% = and J2h-.(h,j)ec x *hj = °- Then > vfc e ^, constraints <\6§ 
can be simplified as follows: 

If < lt k + wn and Zf > Z* fc + w^ - M. 

Hence, the simplified constraints of (J6]) are identical to those of (|2~TI) for Case 1. 
For Case 2, constraints (|2"Tj) can be simplified as follows: 

Zj < l* + w^ - e and 1° > /• + w tj - M. 

According to (!38|) . V/c G P s ,x*^ = and 3fc; G V s , Y^h-(h j)ec x *hj = 1- Then, constraints 
([6]) can be simplified as follows: 

/**' < l'" 1 + w tj - e and if 1 > l- k ' + Wij - M, 

k 



and V/c G V s such that J2h-.(h,j)eC x *hj = 0' 



Zf < l- h + w^ and Zf > l- k + w^ - M. 



Hence, the simplified constraints of ([6]) are identical to those of ( j2~Ti) for Case 2. 
For Case 3, constraints (J2"T1) can be simplified as follows: 

lj < It + Wij and lj > I- + w^. 

According to f[3"8j) . 3ki G P s ,x*-- = 1 and ^h-(h j)^c x * ilj = ■'■■ ^ n addition, according to the 
sub-path optimality constraints, VA; G P s , if cc*^- = 0, X^i:(7i fjec^ftj = 0- Then, constraints 
(jHJ) can be simplified as follows: 

lj h < li l +Wij and lj kl > l kl +Wy, 

and \fk G V s such that x**j- = 0, 

Zf < /• fc + Wy and Zf > lt k + w^ - M. 

Hence, the simplified constraints of ([6]) are identical to those of (1211) for Case 3. 

Since for all the three cases, x*^j,k G V,(i,j) G C results in the same path length 
constraints for DBM as those resulting from y*^, s G S, (i,j) G C for OBM, then, there is a 
corresponding feasible solution satisfying the path length constraints ([6]) in DBM, provided 
that there is a feasible solution satisfying the path length constraints fl2T|) in OBM. D 
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Corollary 2 The path length constraints in OBM are equivalent to those in DBM, concern- 
ing the feasibility of the unique shortest path routing problem. 

Theorem 2 OBM and DBM and equivalent concerning the feasibility of the unique shortest 
path routing problem. 

Proof. By Proposition [TJ DBM is equivalent to DBM'. 

In addition, as discussed at the begin of Section HJ OBM is a reduction of ROBM and 
DBM' is a reduction of RDBM. Besides the additional link weight variables and path length 
variables, the difference between OBM and ROBM are the path length constraints (|2"T|) and 
the difference between DBM' and RDBM are the path length constraints ([6]). 

By Theorem [TJ ROBM and RDBM are equivalent concerning the feasibility of the re- 
laxed problem. By Corollary [2j the path length constraints in OBM are equivalent to the 
counterparts in DBM, and so those in DBM'. Therefore, OBM is equivalent to DBM', and 
so DBM, concerning the feasibility of the unique shortest path routing problem. □ 

Theorem 3 OBM and DBM and equivalent concerning the optimality of the unique shortest 
path routing problem. 

Proof. The conclusion follows directly from Theorem |2]by constructing the corresponding 
optimal solutions between OBM and DBM. D 

Corollary 3 OBM is a correct model of the unique shortest path routing problem. 

5. Comparisons between the Two Formulations 

Concerning the unique shortest path routing problem, it has been shown that routing per- 
formances resulting from the proposed complete formulations are much better than those 
derived from the default methods, by testing on 30 randomly generated data instances with 
combinations of different parameter scenarios. The resulting average maximum utilization is 
30.94% of that from using the hop-count method and 45.54% of that from using the inv-cap 
method. It hence demonstrates the significant gain achieved by formulating the problem 
completely and solving it optimally. 

Between the two complete formulations, compared with DBM, OBM has advantages on 
both constraint structure for applying constraint generation algorithms and model size. 
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5.1 Constraint Structure 

The constraint structures of DBM and OBM are shown in Figure[3]and Figured! respectively. 
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Figure 3: Constraint structure of DBM 

In Figure [31 the first row represents the link capacity constraints §5§, the next four rows 
correspond to the flow conservation constraints (J3J), and the last four rows represent the path 
length constraints (l6l). Accordingly, columns correspond to variables. 

As can be seen, among the three sets of constraints, the flow conservation constraints 
and the link capacity constraints contain only the routing decision variables, whereas the 
path length constraints couple the routing decision variables with the link weight variables 
and the path length variables. Hence, constraint generation algorithms such as the Benders 
decomposition method (Benders, 1962) may be considered to be the promising solution 
approaches for the problem. The problem can be decomposed into one integer programming 
master problem and one linear programming subproblem. The master problem deals with the 
flow conservation constraints and the link capacity constraints, and so contains the routing 
decision variables only. Accordingly, the subproblem copes with the path length constraints. 
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Figure 4: Constraint structure of OBM 



Similarly, in Figure HJ the first four rows represent the path uniqueness constraints ( fl9l) , 
the next four rows correspond to the flow bound constraints ( TT7j) . the third four rows rep- 
resent the flow conservation constraints ffl6l) . the next row corresponds to the link capacity 
constraints ( !i~8l) . and the last four rows represent the path length constraints (121]) . Columns 
correspond to variables accordingly. 

As can be noted, although DBM has a simpler constraint structure, OBM has more 
flexibility to apply decomposition algorithms to solve the problem. 

As shown in Figure HI with OBM, the problem can be globally decomposed into one 
master problem and two subproblems, instead of one master problem and one subproblem 
as with DBM. The master problem contains only the routing decision variables and the 
path uniqueness constraints accordingly. The first subproblem deals with the auxiliary flow 
variables and the second subproblem copes with the link weight variables and the path length 
variables. In addition, the master problem can be further decomposed, with one independent 
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subproblem corresponding to each origin node. 

5.2 Model Size 

Compared with DBM, OBM defines explicitly the auxiliary flow variables and the flow 
bound constraints accordingly. However, in general, |«S| << \T>\, and the size of OBM is 
much smaller than that of DBM. The model sizes of the two formulations are as shown in 
Tabled], where # Variables represents the number of variables and ^Constraints denotes the 
number of constraints. 

More concretely, the model sizes of both the original problems and the master problems 
of the two formulations on a randomly generated data instance with \J\f\ = 50, \C\ = 642, 
\T>\ = 1000, and |«S| = 50 are shown in Tabled 

As can be seen from Table [2J with OBM, the number of variables of the original prob- 
lem decreases from over 600, 000 to 64, 200 and the number of constraints drops from over 
1, 000, 000 to less than 38, 000. In addition, with OBM, both the number of variables and 
the number of constraints of the master problem decline 20 times. 

Table 1: Model sizes of DBM and OBM 

Model ^Variables ^Constraints 

DBM \V\\£\ + \S\\N\ + \£\ |P||jV| + 2|X>||£| + |£| 

OBM 2|S| £| + |S| JV| + £| 2|S||JV| + 3|«S||£| + |£| 



Table 2: Model sizes of DBM and OBM of a large data instance 


Original Problem 


Master Problem 


^Variables ^Constraints 


^Variables ^Constraints 


DBM 645,142 1,334,642 
OBM 64, 200 37, 742 


642, 000 50, 642 
32,100 2,500 



As a conclusion, compared with DBM, OBM has a smaller model size and a more flexible 
constraint structure for decomposition algorithms such as the Benders decomposition method 
to solve the problem. 

6. Conclusions 

With the aim of an exact solution approach to the unique shortest path routing problem 
on average data instances arising from real-world applications, two complete and explicit 
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mathematical formulations with a polynomial number of constraints for the problem are 
developed. A demand-based formulation is first introduced, based on the study of the re- 
lationships between the length of a shortest path and the weights of links that the path 
traverses. The problem is further formulated as an origin-based model by analyzing solution 
properties of the problem. The two formulations are then mathematically proved to be cor- 
rect and to be equivalent concerning both the feasibility and the optimality of the problem. 
Based on the study of the constraint structures and model sizes of the two formulations, the 
origin-based formulation is identified to be the better one for decomposition algorithms such 
as the Benders decomposition method to solve the problem. 

The two formulations may be generalized to other network flow and network routing 
problems. Prospective future work may lie in investigating possible improvements concerning 
both problem formulation and solution algorithm to improve the efficiency of the solution 
approach proposed. In particular, investigations may focus on possible improvements related 
to three factors: the closeness between the initial solution and the final solution to the master 
problem, the strength of cuts generated at each iteration, and the efficiency of an algorithm 
to solve the integer programming master problem. For example, redundant constraints may 
be generated to tighten the feasible region of the initial master problem, strategies such as 
active set method may be applied to strengthen the cuts generated from the subproblems, 
and schemes such as the Lagrangian relaxation method may be embedded into the solution 
algorithm to improve the efficiency of solving the master problem at each iteration. 
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